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Abstract 

We discuss coherent and incoherent (j) meson photoproduction off the deuteron at low energy and 
smah momentum transfer with the aim to check whether the recent experimental data need for 
their interpretation the inclusion of exotic channels. Our analysis of the differential cross section 
and spin-density matrix elements shows that the existing data may be understood on the base of 
conventional dynamics. For a firm conclusion about a possible manifestation of exotic channels 
one has to improve the resolution of the data with providing additional information on channels 
with spin- and double-spin flip transitions being sensitive to the properties of the photoproduction 
amplitude in 7p and reactions. 
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I. INTRODUCTION 



The investigation of the 0-meson photoproduction at low energies, ~ 1.6 — 3 GeV, 
plays an important role in understanding the non-perturbative Pomeron exchange dynam- 
ics and the nature of (t)N interaction. It was expected that in the diffractive region the 
dominant contribution comes from the Pomeron exchange, since the processes associated 
with conventional meson (quark) exchanges are suppressed by the OZI rule [1~7]. An ex- 
ample of such a (suppressed) process is the pseudoscalar tt and 77 meson exchange which, as 
a rule, were considered as a small correction to the dominant Pomeron exchange channel. 
The Pomeron exchange amplitude is usually described in terms of the Donnachie-Landshoff 
model [8], where the Pomeron couples to single constituent quarks as a C = isoscalar 
photon or/and its two-gluon exchange modification [6, 9, 10]. These models are designed for 
the vector meson photoproduction at high energy and small momentum transfer. The valid- 
ity of an extrapolation of these models into the low energy region and close to the threshold 
is not clear. Near threshold, the models predict a monotonic increase of the differential 
cross section of 7p — > 0p reaction at forward photoproduction angle with energy. However, 
a recent analysis of the photoproduction at low energy by the LEPS collaboration shows 
a sizeable deviation from this prediction, in particular, the data show a bump structure 
around E"^ ^ 2 GcV [11]. Another peculiarity of the LEPS data is a strong deviation of 
the spin-density matrix element p\_i from 0.5, which is in favor of a sizable contribution of 
un-natural parity exchange processes. These facts rise several questions: (i) whether one has 
to modify the conventional Pomeron exchange model at low energy, (ii) what is the source of 
un-natural parity exchange channels, (iii) whether we need to introduce some exotic chan- 
nels (additional Regge trajectories, processes associated with possible hidden strangeness in 
the nucleon, etc.) to describe the data. In principle, these questions are related to each 
other and have to be analyzed simultaneously. Thus, for example, the mentioned bump- 
like behavior may be a result of the interplay of the pseudoscalar exchange amplitude and 
modified Pomeron exchange channels. 

The coherent photoproduction off the deuteron in the diffraction region seems to be 
very useful for such an analysis. First of all, the isovector 7r-meson exchange amplitude 
is eliminated in case of the isoscalar target. Therefore, the appearance of the bump-like 
structure in the energy dependence of the differential cross section of the reaction — > (j)D 
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would favor a modification of the conventional Pomeron exchange amplitude. The next step 
is an analysis of spin observables, in particular, the properties of the decay — > K'^K~ with 
unpolarized and polarized photon beams. The incoherent photoproduction in (ppn 
reaction allows to extract observables of the reaction 'jn (pn which can be used for a 
simultaneous analysis of photoproduction off neutron and proton targets in order to get 
additional and independent hint to a manifestation of possible exotic channels. 

Schematically, the coherent and in-coherent (p meson photoproduction processes are ex- 
hibited in Fig. 1 (a,b) and (c,d), respectively. The single and double scattering diagrams 




(c) (d) 

FIG. 1: Diagrammatic representation of coherent (a,b) and in-coherent (c,d) (f) meson photopro- 
duction in 7D reactions with single (a,c) and double (b,d) scattering contributions. 

are shown in (a,b) and (c,d), respectively. The internal dashed line in (b) and (d) corre- 
sponds to "diagonal" (m = 0) and "non-diagonal" (m = TT,p,uj ...) transitions, respectively. 
In this paper wc study the (p meson photoproduction at low energies with < 3 GeV at 
forward photoproduction angles with momentum transfer \t\ < 0.4 GeV^, where the single 
scattering processes are dominant. The coherent 0- meson photoproduction at higher values 
of |i| is controlled by the double scattering processes, which can provide important informa- 
tion about the cross section of the cpN scattering [12, 13]. However, this interesting topic 
is beyond scope of our present analysis, where we focus just on the extremely forward (p 
meson photoproduction, where some hint to an "anomaly" in the differential cross section 
of 7p — > (pp reaction was found [11]. Some theoretical estimate for the coherent vector 
meson photoproduction from deuteron is given in Ref . [14] . The first experimental data on 
— > (pD reaction are reported recently in Refs. [15, 16]. 

The aim of the present paper is to extend the results of Ref. [14] for the coherent and 

incoherent (p meson photoproduction off the deuteron and give a consistent analysis of the 
recent experimental data towards understanding whether they can be described in terms of 
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conventional dynamics or one needs to introduce some new (exotic) processes. 

Our paper is organized as follows. In Sec. II we provide equations for the amplitudes of 
photoproduction off the proton which are used later on for coherent and incoherent (f) meson 
photoproduction in 7Z) reactions. Here we also analyze the unpolarizcd differential cross 
section of the reaction 7^ (pp. In Sec. Ill we present a model of the coherent 'fD — > (pD 
reaction. The incoherent jD (pnp reaction is considered in Sec. IV. In Sec. V we provide 
a simultaneous analysis of spin-density matrix elements for (p — > K^K~ decay distributions 
in 7p, 7n, and reactions. The summary is given in Sec. VI. 

II. * MESON PHOTOPRODUCTION OFF THE PROTON 

For the reaction 7p (pp, we define the kinematical variables with usual notation. The 
four-momenta of the incoming photon, outgoing vector meson, initial and final protons are 
denoted &s k-y, q^, p and p', respectively. The standard Mandelstam variables are defined as 
t^{p'- pf = {k, - q^r, s^W'^(p + k,f. 

In forward-angle photoproduction the s and u channels with an intermediate nucleon and 
nucleon resonances are negligibly weak and the main contribution comes from the Pomeron 
and pseudoscalar (tt, t]) meson exchange processes. The corresponding model for the (p 
meson photoproduction in 7p — > (pp reaction is described in Ref. [7]. However, for the sake 
of completeness in this section we provide the main expressions for the invariant amplitudes 
which will be used below. 

The photoproduction amplitude is expressed in standard form 

where and Sx^ are the polarization vectors of the photon and (p meson, respectively, and 
Ui—u^.{p) [uf—Umf(p')] is the Dirac spinor of the nucleon with momentum p \p'] and spin 
projection [irif]. 

For the Pomeron exchange amplitude we utilize the modified Donnachie-Landshoff (DL) 
model [8] to write 

M^" = M(s, t) , (2) 
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where the transition operator F^^'^ reads 

•y/^^l' 7 -1/7 ^ 



with p — {p + p')/2. The last term with p is added to restore the gauge invariance [7]. The 
scalar function Mp{s,t) is described by the Regge parametrization, 

a pit) 



1 / S \ 

Mp{s, t) = Cp F^{t) F^{t) - — exp 

S \ Sp I 



-—ap{t) 



(4) 



where F-^it) is the isoscalar form factor of the nucleon and F2{t) is the form factor for the (j) 
meson-photon-Pomeron coupling [8] 

^ AMI - O'lt ^ 2f4 

'^'^ {AMI t/toY ' (1 - t/Ml){2i,l + M| - t) ■ 

The Pomeron trajectory is known to be ap{t) — 1.08 + 0.25 1. The strength factor Cp is 
given by 

CP = ^. (6) 

where 7^ ~ 6.7 is the meson decay constant. The parameter g"^ is a product of two 
dimensionless coupling constants — gpss ■ gpqq — {-y/spPs) ■ {\/spPu), where gpss and 
gpgq have a meaning of the Pomeron coupling with the strange quark in (f) meson and light 
quark in a proton, respectively. In our study we choose: to — 0.7 GeV^, — 1.1 GeV^, 
Sp = A GeV^, Ps = l-AA and (3u{d) =2.04 GeV~^ The parameter oat = 2 is taken to be 
larger than the corresponding parameter in DL model [8], making the overall form factor 
close to that of the two-gluon exchange model [10]. Actually, the original DL model was 
motivated by the two-gluon exchange model of Landshoff and Nachtmann [17], therefore 
such a modification seems to be reasonable. 

In the case of the pseudoscalar mesons exchange (M = tt, 77), the transition operator Al^^ 
reads 

M,^ = -^ j^^ 75 ^_M2 ^mW (7) 

with gT^NN — 13.26, gi-y,^^ ~ —0.14, and gi-y^^ ~ —0.71 [7]. In this paper, following estimates 
based on QCD sum rule [18] and chiral perturbation theory [19], as well as the phenomeno- 
logical analysis of rj photoproduction [21], we use gr,NN — 1-94. is the product of the 



two form factors of the virtual exchanged mesons in the MNN and ^VM vertices 

FMit) = (8) 



with A^(^) = 1.05 GeV. This value is shghtly greater than the values of cut-off parameters 
in Ref. [7] (A^(^) = 0.6 (0.9) GeV), which result in some modification of the pseudoscalar 
exchange contribution. The SU(3) symmetry predicts a constructive tt — rj interference in 
7P reactions and a destructive interference in ^yn reactions [22] . 
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FIG. 2: Differential cross section of the (j)p reaction as a function of momentum transfer t 

at = 2.02 GeV. The Pomcron, pseudoscalar exchange contributions and the total cross section 
are shown by dot-dashed, dashed and solid curves, respectively. Circles, triangles down, triangles 
up and squares correspond to the LEPS [11], SAPHIR [25], Bonn [26], and JLab [27] data, 
respectively. 

In Fig. 2 we show the differential cross section of the 7p — > (j)p reaction (sohd curve) for 
the photon energy bin Ej — 1.97 — 2.07 GeV from LEPS [11], together with the experimental 
data at Ey 2 GeV [11, 25, 26]. For completeness, we also display JLab [27] data, obtained 
at 3.6 GcV, because there is no much difference in the t dependence of Bonn [26] and 
JLab [27] data. One can see that the model satisfactorily describes the Bonn and JLab 
experimental data. However, it underestimated the LEPS and SAPHIR data at relatively 
large \t\ which probably may manifest additional channels beyond our simple model [7]. The 
energy dependence of the differential cross section at forward photoproduction angle with 
^ = (i.e. t — tmax) together with the experimental data [11] is shown in Fig. 3. One can 
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see a sizeable deviation of experimental data around — 2 — 2.3 GeV from the monotonia 
theoretical curve, which is related to the difference between ^-dependence of different data 
sets and our model, discussed above. It is clear that for understanding the nature of this 

difference one needs more precise experimental data not only in differential cross section but 
in polarization observables sensitive to the spin flip channels at several energies. 

III. * MESON PHOTOPRODUCTION IN REACTIONS 

In this section we consider coherent 7D — > (pD and incoherent 7D — > (pnp photoproduc- 
tion processes. The kinematical variables for these reactions are the following ones. The 
four-momenta of the initial and the final deuteron {np system) are denoted as po and p'x 

[X = D,np), respectively. The Mandelstam variables are defined as sd = W^, = {po + k^)^, 
= {p'x ~ PdY 1 and so on. The space component of the momentum transfer to deuteron 
in the laboratory system is = — —toi'^ — to/^My), where Md is the deuteron mass. 

A. Coherent photoproduction 

As mentioned above, here we consider the meson photoproduction at forward angles 
with \t\ < 0.4 GeV^, where the dominant contribution comes from the single scattering 
process, shown in Fig. 1 a. In such a case one can use a non-rclativistic framework for the 
deuteron form factor based on utilizing the realistic NN interaction. In our analysis we 
use the deuteron wave function calculated with Paris potential [23, 24] designed just for 
describing nuclear processes at high momentum transfer. Thus, it describes fairly well the 
deuteron electromagnetic form factor with momentum transfer up to —t ~ 0.9 GeV^ [24]. 
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FIG. 3: Differential cross section of the 7p — > 



reaction at t 



= imax {0 = 0) as a function 



of the photon energy. The experimental data 
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are taken from [11, 28]. 
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The total vector meson photoproduction amplitude in the reaction 7D — > VD reads 

a/3 

where Mj, Mf, A-y, and Ay stand for the deuteron-spin projections of the initial and final 
states, and helicities of the incoming photon and the outgoing vector meson, respectively, 
is the amplitude of the vector meson photoproduction from the isoscalar nucleon 

r^ = -(TP + r"). (10) 

The indices a and P in Eq. (9) refer to all quantum numbers before and after the collision. 
The "elementary" photoproduction amplitudes T^'" are defined in the previous section, tt 
exchange terms are canceled in the total amphtude since = 

Using the standard decomposition of the deuteron state in terms of s (Uo) and d {U2) 
wave functions, one can rewrite Eq. (9) in the explicit form 

_ o^fiz\^ a\^'^y. (n^r^^ r^^"^' r^^^^ r^^^f 

^ C'l'M^, A//^L'OAO RLL'\{q^) ^mim'i;AvA^ (0 ' (H) 

where j — ^J2j + 1, and the radial integral Rll'x reads 

RLL'x{q^) = j drUL{r)UL'{r)jx{qr/2). (12) 

For a qualitative analysis of the unpolarized differential cross section at small momentum 
transfer with ~ 0, keeping only the spin/helicity conserving terms with natural and 
unnatural parity exchange in the total amplitude, one gets 

N / 1 \ N 

T^m'-MX^it) = [^^^J Smm'Sx^X^T- (t). (13) 

N 

Here, Tq {t) is the spin-independent part of the amplitudes. Using Eq. (11) with Eq. (13), 
we get the following result for the natural and un-natural parity-exchange parts of the total 
amplitude 

-<DN ox X (a cN I r cN\rnN 



'^MfMi-MX^ — '^^MiMfSx.yXv{^±lMiSi + ^OMi'S'o )Tq , 

^M^M,;AvA^ = '^Mi\y6MiMfSx^Xy5±iMi . (14) 
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The form factors Sf'^ read 



^ Fc- ^/2Fq, ^Fc + 2V2Fq, = Fm , (15) 



with 



Fc — -Rooo + -R22O, Fq — R202 7^-^220 ) 



1 

71^ 



Fm — Rooo — 2^220 + V2R202 + R220 ■ (16) 

Taking into account the cancelation of the un-natural parity tt exchange contribution and 
neglecting weak 77 meson exchange, one can express the differential cross section of the 
(f)D reaction by the cross section of the (j) photoproduction from the isoscalar nucleon 

< > as 

-IF ^ ' ^''^ 

where t — to and Z{t) is the structure factor 

Z{t)^F'a{t)+^F'Q{t). (18) 




FIG. 4: The dependence of the structure fac- 
tor Z on t = tn- 
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The dependence of Z and Fc,q on to is rather symbohc. In fact, these factors depend 
on the spatial part of the four-momentum transfer in the laboratory system q, as follows 
from Eq. (12). The relation between tr, and reads tf) = —2M£)(^y q"^ + M|, — Md). The 
structure factor Z as a function on tD is shown in Fig. 4. In the considered region of 
momentum transfer t, the factor Z{t) is related to the well known structure function A{t) 
of the elastic eD — > eD scattering as 

A{t)cZ{t)Gl{t), (19) 
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where Gd{t) = 1/(1 — i/0.71)^ is the dipole electromagnetic form factor of the proton. 

Equation (17) allows to "extract" the cross section of the ^ < N > reaction from the 
measured cross section of the ^yD reaction as 



Jt 



dt 



(20) 
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FIG. 5: Differential cross section of the 
7Z) (j)D reaction as a function of momen- 
tum transfer t {t = tD)- Circles and squares 
correspond to LEPS [15], and CLAS [16] 
data, respectively. 



In Fig. 5 the differential cross section of (f)D reaction is exhibited calculated by 

using the explicit expression for the photoproduction amplitude given by Eq. (9), together 
with the available experimental data by LEPS (circles [15]) and CLAS (squared [16]) col- 
laborations. For simplicity, we show only a comparison for the bin — 2.07 — 2.17 GeV. 
The description of the data for other bins has a similar quality. One can see that the model 
rather well describe the the data at low momentum transfers \tD\ but tends to underestimate 
the data at higher probably pointing to growing weight of more complicated (such as 
double scattering) channels. 

In Fig. 6 we show the energy dependence of the differential cross section of the — >■ (pD 
reaction at 6* = (i.e. t — imax) together with experimental data [15]. The agreement 
between data and model is fairly reasonable. Note that here the experimental data do not 

2 GeV. 



point to a bump-like structure at 
10^ 




FIG. 6: Differential cross section of the 
7Z) — ^ (l)D reaction at i = tmax {i = to) 
as a function of the photon energy. The ex- 
perimental data are taken from [15]. 



In Fig. 7 the comparison of meson photoproduction off the proton and off the isoscalar 
nucleon in a deuteron at ^ = is displayed. In latter case the experimental data and 
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FIG. 7: Differential cross section of meson 
photoproduction off the proton (solid curve) 

and off the iso-scalar nucleon in 7D reaction. 
The experimental data are taken from [11, 
15]. 

the theoretical curve are evaluated from the corresponding cross section of the — > (pD 
reaction by using Eq. (20). The figure displays the energy dependence of the differential 
cross sections at ^ = 0. One can see that the two cross sections are close to each other 
at all energies. The Pomeron exchange amplitude dominates ah high energies. At lower 
energy, the behavior of cross sections of the 7p and ^ < N > reactions is not trivial. The 
elimination of the isovector tt exchange contribution in the < N > reaction is compensated 
by a modification of momentum transfer t, which is smaller compared to that of 7p reaction 
near the threshold in jD reaction. This causes the approach of both curves with decreasing 
energy E^. 




B. Incoherent photoproduction 

The main purpose of the measurement and the theoretical study of the incoherent 0- 
meson photoproduction in ^yD reactions is an extraction of the cross section of jn 4>n 
photoproduction with the goal of a subsequent combined analysis of 7p and 772 reactions to 
seek for a possible manifestation of exotic channels. This problem seems not too difficult if 
one uses the exclusive jD (pnp reaction. But at low energy and forward photoproduction 
angles, the momenta of the recoil nucleons are small, and there is an experimental problem 
with their detection. Therefore, another way is to study the [jD, (f)] missing mass distribution 
in the inclusive jD 0X (X = np, D) reaction. Below we develop a model which can be 
used for an extraction of the observables of 7n (j)n photoproduction. 

The differential cross section of the (p meson photoproduction in the 7D — > (pnp reaction 
reads 

da 1 f j?^ P /\rn |2 



dtdMx 167r(s -M|,)2y IGtt^ 



d^T^,i\TpV + \Tn\'). (21) 
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where p and Q are the momentum and the sohd angle of the spectator nucleon in the rest 
frame of the np pair, respectively; Mx is the invariant mass of this pair, and t — tx] averaging 
and summing over the spin projections in the initial and the final states are assumed. Tp(„) 
is the amplitude of the partial proton (neutron) contribution. It is related to the amplitude 
of the <pN [N = n, p) reaction and the deuteron wave function il)^ as 



= J](-m2-m|lM, - A)(LA1M, - A|1M,) T^f,;^^ ^£v(p.) (22) 



where Ps is the spectator momentum in the laboratory system, ULij-) is the radial deuteron 
wave function in the configuration space, Mj A^, mi^2) and \^ are the spin projections of the 
incoming deuteron, photon helicity, the spin projections of the outgoing nucleons and the 
helicity of the (j) meson, respectively. For evaluating Eq. (21) we define kinematical variables 
by the following steps. For given Mx, the energy of the outgoing nucleons in the np rest 
frame is £^ = Mx/2. Then, using Q and the 0- meson photoproduction angle in the center 
of mass system as input variables we evaluate the four-momenta of the outgoing nucleons 
first in c.m.s. and then in the laboratory system. The four-momentum of the struck nucleon 
is, Pi = pd — Ps, where Pd = (^d,0). The amplitude T''^ in Eq. (11) is evaluated with 
an off-shell struck nucleon with < < M^^. In such a way, the off-shell effects in the 
incoherent channel are evaluated consistently. 

The differential cross sections of the incoherent meson photoproduction are displayed 
in Figs. 8 and 9. Let us first discuss the differential missing mass distribution in the 
— » (pX reaction (X = D,np) as a function of the ['~fD,(f)] missing mass and momen- 
tum transfer t. For the coherent and incoherent parts we use the common momentum 
transfer t — to- This means that the incoherent part must be multiply by the Jacobian 



dtx/dto — Y A(s, M|-, M^)/X{s, M|,, M|). With regards to a comparison of our prediction 
to the experimental data, the experimental resolution must be included. Also, the cross 
section of the incoherent photoproduction is slightly modified. Therefore, we compare data 
with the missing mass distribution folded with a Gaussian distribution function 



LA 



with 



V^f^(p) = {2n)h'^YL^{p)uL{p) 




(23) 
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FIG. 8: Differential distribution of [yD, cp] missing mass in jD (j)X reactions at different energies. 
The curves correspond to the cross sections of incoherent 7D — (t>{pn) (dashed) reactions, the 
coherent 7D jD (dot-dashed), and their sum (soUd). The arrows mark the position of the 
maximum of the missing mass distribution. 



f{Mx - M) 



exp 
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(24) 



with (T = 10 MeV [15], which imitates a finite experimental resolution. 

In Fig. 8 we show the differential [7D, 0] missing mass distribution in (pX reactions 

for different photon energies for forward photoproduction angle = 0. The position of the 
maximum of the incoherent part is marked by an arrow. One can sec a strong energy 
dependence of (i) the absolute value of the cross section, (ii) the relative contributions of the 
coherent and incoherent processes, (iii) the position of the maximum of the incoherent part. 
At relatively large photon energies {Ej ~ 2.5 GeV) our model predicts a strong overlap of 
coherent and incoherent parts, and the coherent photoproduction amounts more than 30% 
of the total cross section. Our model seems to be an effective tool to isolate the coherent 
and incoherent parts with subsequent extraction of the (p photoproduction off the neutron. 



Fig. 9 exhibits the invariant mass distribution averaged within the interval 



1.5 - 



2.4 GeV together with experimental data [29] given in units of events. The theoretical curves 
are scaled by the factor 3.7 [//b/GeV^]~^. The comparison is rather qualitative because we 
did not use the detailed acceptance corrections which may somehow modify the shape of the 
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distributions. Nevertheless, the quahtative agreement between prediction and data seems 
to be quite encouraging. 



IV. SPIN DENSITY MATRIX ELEMENTS 

In this section we consider several important matrix elements of spin-density matrices 
p\\' ~ 0? 1?2) which determine the meson decay distribution in its rest frame in case 
of both unpolarized and linearly polarized photon beams. The spin-density matrices are 
defined by 



The symbol a includes the polarizations of the incoming and outgoing baryons, and the 
normalization factor has the standard form 



where T^-^x^x-, is the total meson photoproduction amplitude. 

We perform our consideration in the 0-meson rest frame with the quantization axis along 
the beam momentum, i.e. Gottfried- Jackson (GJ) system. Other possible choices are the 
helicity (H) system with quantization axis opposite to the recoil nucleon (dcuteron) momen- 
tum in 7p reaction, and the Adair (A) system, where the quantization axis is along 
the beam direction in c.m.s. [30]. The GJ system has some advantage because only here 
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FIG. 9: Distribution of [7Z), 0] missing mass 
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corresponds to the experimental data [29]. 
The theoretical curves are scaled by the fac- 
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tor 3.7 [//b/GeV^] ^ (see details in text). 
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some of spin-density matrix elements have a clear physical meaning, e.g. as a measure of 
the helicity conserving processes or as an asymmetry between processes with natural and 
un-natural parity exchange in t channel. 

Consider first the matrix element Pqq. This matrix element determines the polar angular 
distribution of — > KK decay 

3 



w{cose) 2 



(poo + ^(l-3p;)sin^e) . (27) 



In GJ system, Pqq is the measure of the spin flip transition with A-y = ±1 — > A<^ = 0. Thus, 
in case of a pure helicity conserving amplitude, which may be expressed as 

Ta;X4„\^ - (£a^ " ^A^) , (28) 

the photon polarization vector e;^^ is transversal with respect to the z axis, and therefore 
spin-flip transitions A-y = ±1 — > = are forbidden and Pqq = 0, independently of the 
momentum transfer. In the helicity system, the photon polarization vector has a finite z 
component 

£A,. = ^sin/3, (29) 

where (3 is the angle between H and GJ systems 

0=^^^, (30) 
COS B — 1 

and and 9 are the 0-meson velocity and the photoproduction angle in c.m.s., respectively. 
Foe relatively large momentum transfer, when sin /? ~ 1, one gets a large value of 

~ sin^ (3 , (31) 

even for the helicity conserving amplitude. Conversely, one can imagine an amplitude which 
generates Poo*^ — 1 (fo^" example, take only the second term in Eq. (3)), and then Pqq ~ 
cos^ ~ 0. In general, the spin-density matrices in H and GJ system are related to each other 
as 

P\''x'-T.dU-P)p',fdUP)- (32) 

Let us first discuss the energy dependence of the spin-density matrix element Pqq in the 
GJ system for 7p, ^yn and 7D reactions. Following the experimental data, we calculate 
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averaged p matrices in the interval \t\ — |io| < ^t- The averaged p matrices are defined as 
ratios of averaged numerators and denominators {N) in Eqs. (25). In such direct 
comparison of the p matrices for the coherent and for the 7p reactions is hampered by 
the deuteron form factor. The deuteron form factor drops rapidly with increasing values 
—t (see Fig. 4) and, therefore, the dominant contributions in and 7p reaction at the 
same values of and Aj come from different momentum transfers I^dI < \ip\. This effect 
is particularly important for small values of ~ ^max, where the slope of the deuteron 
form factor is rather steep. Thus, at relatively large energies, say > 2 GeV, the main 
contribution comes from ~ |to| ~ 0, making the averaged p matrices for the reaction 
practically constant. One can remove the effect of the deuteron form factor by scaling the 
product TT"^ in Eqs. (25) (or in the cross sections of the K^K^D reactions) by an 

inverse structure factor Z{t) given by Eq. (19). Such reduced p matrices would be much 
closer to the p matrices for the photoproduction off the "free" isoscalar nucleon. Fig. 10 
illustrates effect of the deuteron form factor for the case of the reactions and the 7D 




(GeV) (GeV) 

FIG. 10: The energy dependence of Pqq for the -yD (pD reaction, (a) and (b) correspond to 
to = tmax and to = —0.2 GeV^, respectively, = 0.2 GeV^. The solid, dot-dashed and dashed 
curves correspond to the case of the expUcit jD reactions, the 7D reaction with reduced cross 
sections, and photoproduction off the free isoscalar nucleon, respectively 

reaction with reduced cross sections. The latter one is denoted as 7-Dr- For completeness, we 
also show results for the cj) photoproduction off the free isoscalar nucleon. One can see a large 
difference between predictions for reaction and the photoproduction off the free isoscalar 
nucleon at to = tmax- In the first case, Pqq is almost constant, whereas in the second case 
it increases with energy in the given energy interval. Such an increase for the jN reaction 
can be understood as follows. The finite value of Pqq is generated by the Pomeron exchange 
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amplitude and is determined by the second (main) and third terms in Eq. (3), whereas the 
total cross section is dominated by the first term. Neglecting spin conserving pseudoscalar 
meson exchange one can get the following analytical estimate of Pqq for GJ frame for the 
pure Pomeron exchange channel 



approx 



M%){s 



Ml-Ml-t) ' 



(33) 



where Px is the x component of the nucleon momentum {p^ — p^), k-y is the photon energy, 
and s is the total energy squared in the 7A'" vertex. At fixed dependence on form factors in 
numerator and denominator for the reaction is canceled. The increase of Pqq with energy, 
within the considered energy interval, is explained by a faster increase of the numerator 
(because of factor p^) compared to the denominator at fixed t. At larger energies and small 
\t\ this ratio and the corresponding matrix element decrease. 
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FIG. 11: Estimates of Pqq given by 
Eq. (33) for and 7Z) reactions. 



The difference between reduced Pqq matrix element and the case of photoproduction off 
the isoscalar nucleon is explained by the difference in p^, k^, tmax, and s for 77? and ^yD 
reactions. Actually, the kincmatical variables in 7A^ vertices in ^p and 7-D reactions at fixed 

and t (Itf^axl < l^maxl) different and this difference is reflected in spin-density matrix 
elements. As an illustration, in Fig. 11 we exhibit results for given as a ratio of the 
averaged numerator and denominator in Eq. (33) calculated for 7p and jD kinematics. One 
can see some difference between the two cases caused by pure kinematics. 

The comparison of Pqq for the jp, jn and jD reactions without and with scaling by Z~^{t) 
is shown in Fig. 12. In Fig. 12 (a) we show the result for forward photoproduction angles with 
^0 = ^max (^ = 0), together with available experimental data [11]. In Fig. 12 (b) we choose 
the case of a larger momentum transfer with to — —0.2 GeV^ for each energy. One can see 
a monotonic increase of Poo with energy, and the inequality Poo(7P) < Pooil''^) < Pooil-^r) 
holds. Some enhancement of Pqq in 771 reactions is explained by the destructive interference 
in the tt — 77 meson exchange amplitude which leads to a decrease of the helicity conserving 
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terms in the full amplitude. Therefore, the relative contribution of the spin-flip terms in the 
7n reaction (cf. Eq. (3)) would be larger. In case of the ^D,. reaction, together with a total 
suppression of tt meson exchange, Pqq increases additionally because of some difference in 

kinematics, as discussed above. 

In Fig. 13 we exhibit the angular distribution W^(cosO) in the (pD K^K~D 



reaction in the helicity frame for — 3.1 GeV and for to — —0.3 GeV^ together with 
available experimental data [16] given in this frame. The shown experimental data are 
obtained in two energy bins with — 1.6 — 2.6 and 2.6 — 3.6 (GeV)and momentum transfer 
|t| = 0.35 — 0.8 GeV^. In our calculation the momentum transfer is in the range \t\ = 
0.3 — 0.5 GeV^, which corresponds to an upper bound of the momentum transfer acceptable 
for our model for the 'jD — >■ (pD reaction with single scattering processes. Nevertheless, one 
can see a reasonable agreement between calculation and data. Note that this distribution 
is different in different frames because of the frame dependence of the p matrices. As an 
example, in Fig. 14 we show the energy dependence of Pqq for the — > cpD reaction in H 
and GJ frames at \t\ - lU < 0.2 GeV^ and -to = 0.2 GeV^. 
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FIG. 12: The energy dependence of Pqq. (a) and (b) correspond to to = ^max and = —0.2 GeV^, 
respectively. The experimental data are taken from [11]. 
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FIG. 13: The angular distribution M^(cos9) for 
the — ^ (f)D — ^ K'^K'D reaction in the heUc- 
ity frame at i?^ = 3.1 GcV and —Iq = 0.3 GeV^. 
The experimental data for two energy intervals 
and |i| = 0.35 - 0.8 GeV^ are taken from [16]. 
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The energy dependence of the spin-density matrix element Rep^-i is displayed in Fig. 15. 
This matrix element determines the azimuthal angle distribution of — > KK decay in 
reactions with an unpolarized photon beam 



W\^) = —(1 - 2Rep?_iCos2$) . 



(34) 



The matrix element p^^i is proportional to the relative contribution of processes with double 
spin transition where = ±1 — = =Fl. In our model, these transitions are generated 
by the last term in Eq. (3). In Fig. 15 a (b) we show results for \t\ — \to\ < 0.2 GeV^ with 
^0 = ^max (—0.2 GeV^), together with available experimental data [11]. The reason of the 
inequality Pi_i(7p) < Pi-iiin) < p5-i(7-D) is similar to that in the previous case of single 
spin-flip transitions. 

The matrix elements p\'^i are related to the asymmetry of transitions with natural (flrst 
term of Eq. (3)) and un-natural (tt,?]) parity exchange. They determine the meson decay 
distribution in case of linearly polarized photons as a function of the angle between azimuthal 
decay angle ($) and the angle of the polarization plane {^) 



H/^($-^) = — (l + 2P^p}_iCos2($-^)) , 



where is the strength of polarization and 



(35) 



(36) 
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FIG. 14: Spin-density matrix elements pgo 
the helicity and Gottfried- Jackson frames at \t\ — 
\to\ < 0.2 GeV^ and to = -0.2 GeV^. 
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FIG. 15: The same as in Fig. 12, but for Re/9?_i. 
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The energy dependence of the spin-density matrix element p\_i is shown in Fig. 16 
together with the experimental data [11, 15]. In this case, the effect of the deuteron form 
factor is rather weak and we do not display results for the reduced matrix element. For pure 
natural (un-natural) parity exchange it is equal 0.5 (-0.5). Qualitatively, within experimental 
accuracy, the result of our calculation is consistent with the data. Sizeable deviations of p\_i 
from 0.5 in the 7p reaction at low energy is explained by a large contribution of the tt, t) 
exchange processes. Thus, at ~ 2 GeV they contribute on the level of 30% to the 
total cross section. In 771 and ^yD reactions the pseudoscalar exchange contributions are 
suppressed, shifting p\_i towards 0.5. 

For completeness, wc also present the angular distribution W^{^ — \E') of Eq. (35) for 
different cases. Figure 17 exhibits this angular distribution for the reaction 7p ^ 0p — >■ 
pK+K- at 1^1 - Itmaxl < 0.2 GeV^ in two energy intervals = 1.97 - 2.17 and 2.17 - 2.37 
(GeV) with beam polarization — 0.86 and 0.90, respectively, together with available 
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FIG. 16: The same as in Fig. 12, but for p\_i. The experimental data are taken from [11, 15]. 
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FIG. 17: The angular distribution VF^($ — ^') for the reaction ^ <j)p ^ pK^K at |t| — |tmax| < 
0.2 GeV^. (a) and (b) correspond to the energy intervals = 1.97 - 2.17 and 2.17 - 2.37 (GeV), 
respectively. The experimental data are from [31]. 
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experimental data [31]. One can see a reasonable agreement between our calculation and 
the experiment. 

M>,<1.89GeV Mx>1.89GeV 
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FIG. 18: The angular distribution for the reaction jD ^ (f)X ^ XK+K' {X = D, np) 

at \t\ — |tmax| < 0.1 GeV^. (a) and (b) correspond to the [7!^,^] missing mass smaller or larger 
than 1.89 GeV, respectively. The experimental data are from [15]. 

The angular distribution W^{^ — \E') for the inclusive 'jD (j)X {X = D,np) reaction 
is displayed in Fig. 18 together with the experimental data of Ref. [15]. This distribution is 
calculated using the model, developed in Sec. III. The left (a) and right (b) panels correspond 
to events with [7D, 0] missing mass smaller or larger than Mcut — 1-89 GeV, respectively. 
In the first case the contributions come both from the coherent and incoherent (p meson 
photoproduction. The "effective" p\_i matrix element is expressed as a sum 

pl-ies = Pi-iD PcH + Pi-mp (1 - Pch) , (37) 

where Pd is the relative weight of the coherent channel, and Pnp — {pn+pp)/'^ is the p matrix 
for the quasi-free nucleon. In the second case, the contribution of the coherent channel is 
negligible and we get 

pWeS ^ p\-Xnv (38) 

In Fig. 18 we show result for \t\ — |imax| < 0.1 GeV^ and the energy bin with — 2.21 — 
2.37 GeV [15]. Here, the beam polarization is P^ — 0.935 and the model predicts Pch — 0.67. 
One can see a sufficient agreement between the theoretical curves and the data. A similar 
agreement holds for the other energy bins, too. 

The agreement between the experimental data and the calculations for the K^K^ an- 
gular distributions in 7p and ^np reactions means that the model describes correctly the (j) 
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photoproduction off the neutron, and in particularly, supports our choice of the pseudoscalar 
channel with a small contribution of the rj meson exchange. 
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FIG. 19: The energy dependence of 2pli + pJo- The left panel and the right panels correspond to 
*o = ^max and to = —0.2 GeV^, respectively. 

The sum p\j = 2p\^ + Pqq determines the meson decay distribution as a function of the 
angle between production and beam polarization planes 



W^^W) = 2^(1 + 2P^pLcos2*). 



(39) 



It is important that parity conservation requires p^^ — (— [30], which makes 
p]^ invariant under rotation of the coordinate frame in the production plane. This means 
that pj^^ = p™ — p]^. Therefore, it is natural that this invariant function determines the 
distribution which depends only on the beam polarization. 

Since p\j is proportional to a combination of single and double spin-flip transition am- 
plitudes, its absolute value is small. The energy dependence of p\j is shown in Fig. 19. One 
can see some increase of p\j when going from to — tmax (a) to to — —0.2 GeV^ (b). This is 
explained by an increasing contribution of spin- flip transitions with \t\. 

V. SUMMARY AND DISCUSSION 

Wc studied different aspects of coherent and incoherent cf) meson photoproduction off 
the deuteron at forward photoproduction angles with the aim to check whether the recent 
experimental data require the inclusion of some exotic channels discussed in literature. For 
this purpose we re-analyzed the elementary 7p — > (f)p reaction in order to use it as an input 
for our study. The corresponding amplitude in the diffractive region is expressed as a sum of 
Pomeron and pseudoscalar exchange channels. The flrst one represents a shghtly modifled 



22 



Donnachie-Landshoff Pomeron exchange amplitude, whereas the second one is the coherent 
sum of the tt and rj meson exchange channels. In present work the contribution of the rj 
exchange channel is relatively weak, and correspondingly, the tt exchange is enhanced in 
order to get the proper relative contributions of the channels with natural and un-natural 
parity exchange. The Donnachie-Landshoff model is designed for high energy and it is not 
clear whether it can be applied at low energies, and close to the threshold as well. 

We performed a detailed analysis of the differential cross section of the 7p (f)p reaction 
at £^-y ~ 2 GeV and obtained a reasonable agreement between the model predictions and 
the available experimental data in diffraction region. At larger momentum transfer our 
model underestimates recent data of LEPS and SAPHIR but is quite reasonable for the 
Bonn and JLab data up to t = 0.8 GeV^. On the other hand, the Pomeron exchange 
model, motivated by the two-gluon dynamics, contains terms responsible for single and 
double spin-flip transitions. The model predictions for the spin-density matrix elements 
being sensitive to the spin-flip transitions are in agreement with available data for the 7p 
reaction at ~ 2 GeV, which also decreases the space left for possible exotic channels. 
Therefore we can conclude, that for a clear understanding a possible manifestation of an 
exotic channel one needs a complete set of t dependences for unpolarized cross sections and 
polarization observables at different energies. 

We developed a model for the coherent and incoherent meson photoproduction off 
the deuteron and performed again a detailed analysis of the existing data. The slope of 
the differential cross section of the coherent meson photoproduction is deflned by the 
corresponding slope of the elementary 7A'^ reaction and by the deuteron form factor. We 
found a quite reasonable agreement between the model prediction and the experimental data 
in the diffractive region and some underestimate at large |t| ~ 0.4 GeV^, which favor the 
contributions of more complicated channels, for example, double scattering processes. But 
on the other hand, the model calculation of the (f) — > K'^K~ decay distribution, l^(cos©), 
at \t\ ~ 0.4 GeV is in a good agreement with the experimental data, which, to some extent, 
support the single scattering model in this region of t. Therefore, the remaining difference 
between theory and experiment at \t\ ~ 0.4 GcV^ requires further investigation. 

The model fairly well describes the energy dependence of the cross section of the 
(pD reaction at ^ = without any hint to a bump-like behavior. 

We performed detailed and combined investigation of several important spin density 
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matrix elements for 7p — > 0p, coherent 7D — > and incoherent 7D — > 0np reactions 
aimed at (i) studying effect of ehmination of the isovector tt meson exchange in coherent 
reaction, and (ii) extracting observables for the 771 reaction. The ehmination of the tt 
meson exchange has two consequences. One is the relative decrease of channels with spin- 
conserving amplitudes, which result in an increase of the relative contributions of the spin 
flip transitions. This leads to an enhancement of the corresponding spin density matrix 
elements. Another one is related to a strong suppression of the amplitude with un-natural 
parity exchange and shift p\_-^ matrix element towards 0.5. We got a common description 
of (f) meson decay distributions for 7p — > 0p and incoherent — > 0np reactions confirming 
the reliability of our model for the 772 reaction. 

To summarize we can conclude, that the existing experimental data (including also very 
recent data) on 7p, coherent —>■ (pD, and incoherent jD (pnp reactions in the diffrac- 
tion region at low energies support the model based on the dominance of the Donnachie- 
Landshoff Pomeron plus tt, rj exchange channels with a relatively weak 77 meson contribution. 
For a definite conclusion about a possible manifestation of exotic channels one has to im- 
prove the resolution of the data with providing additional information on the channels with 
spin- and double-spin flip transitions being sensitive to properties of the photoproduction 
amplitude in 7p and jD reaction. This problem may be studied experimentally at the elec- 
tron and photon facihties at LEPS of SPring-8, JLab, Crystal-Barrel of ELSA, and GRAAL 
of ESRF. 
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